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Abstract 

We examine the stability of intense flat bunches in barrier buckets used in the Recycler. 
We consider some common stationary distributions and show that they would be unstable 
against rigid dipole oscillations. We then discuss an analytical model for the line density that 
best describes measured bunch profiles. We include space charge in this model to predict the 
bunch intensity at which Landau damping would be lost. The dependence of this threshold on 
the bunch length is studied and related to the results of an experimental study with shorter 
bunch lengths. The threshold for the microwave instability is estimated. These studies will 
be followed by more detailed numerical studies. 
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I. INTRODUCTION 



The creation of long fiat bunches is under study for the LHC upgrade as a way of 
increasing the luminosity The stability of such bunches is one of the key issues of 
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interest. At the Fermilab Recycler {2], long flat bunches are created using rf barriers. 
At present intensities these bunches are observed to be stable with lifetimes around 
20-50 hours (depending on intensity). When electron cooling is enabled, the lifetime 
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Parameter 


Value 


Units 


Circumference 


3331.0 


m 


Energy 


8.93 


GeV 


Bunch intensity 


4.5x10^2 




It 


19.97 




Rf voltage Vq 


1.8 


kV 


Ti 


48.0 


0.0189/Lf sec 


T2 


324.0 


0.0189// sec 


Bunch area 


100.0 


eV-sec 


max 


9.97x10-6 





TABLE I: Recycler bunch parameters 



reaches ~ 500 hours at bunch intensities below 4.5x10^^. We explore the longitudinal 
stability of these bunches at higher intensities [3|, 

Bunches are confined within a barrier bucket by two voltage pulses of equal mag- 
nitude and opposite polarity. The pulses, of equal duration Ti, are separated by a 
duration T2 with no applied voltage. The body of the bunch is contained within the 
interval T2 but the head and tail of the bunch penetrate into the barrier on either 
side. Between the barrier pulses, the beam particles feel no longitudinal focusing force 
and can be considered to be coasting. Changing T2 adiabatically changes the bunch 
length while preserving the bunch area. The main longitudinal parameters of the Re- 
cycler bunches are shown in Table [H Single particle dynamics within such a bucket 



was studied in 



Co 



lective effects were numerically studied in the context of a low 
energy heavy ion ring gJ. Longitudinal stability in the Recycler may be influenced by 
several factors. For example, the synchrotron period is rather long, hence even slowly 
growing instabilities which normally be would not be of concern can in this case be 
important. Furthermore, since the bunches are long, they can be excited by relatively 
low frequency excitations. 

In this paper we will consider primarily the thresholds for the loss of Landau damp- 
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ing and also discuss the results of experiments to determine the impact of a local 
instability possibly induced at short bunch lengths. In Section II we discuss single 
particle dynamics in a barrier bucket for the sake of completeness and handy reference. 
In Section III we consider two stationary distributions that are widely used to describe 
bunches confined by conventional harmonic rf systems as candidates to describe the 
Recycler bunches. In Section IV we calculate the coherent frequency of rigid dipole 
oscillations for these distributions and the relationship to the maximum of the bare 
incoherent frequency. In Section V we find an analytical distribution that adequately 
describes the measured profiles. This distribution is then used to calculate the incoher- 
ent frequency distribution in the presence of space charge. The intensity threshold at 
which Landau damping would be lost is predicted while keeping other bunch param- 
eters unchanged. In Section VI we consider the experimental results obtained when 
the bunch length is changed and discuss the dependence of the threshold for loss of 
Landau damping. Section VII has a brief discussion of the microwave instability and a 
rough estimate of the threshold intensity. We conclude in Section VIII. In a subsequent 
paper we plan to a) discuss the stability diagrams for the Recycler bunches and b) use 
numerical simulations to study the stability limits in greater detail. 

II. SINGLE PARTICLE DYNAMICS 

As mentioned above, a barrier bucket consists of two voltage pulses of equal and 
opposite polarity, each lasting for time Ti and separated by a time T2. The pulses 
can be of arbitrary shape, rectangular, triangular, half-sinusoidal etc. In most of this 
report we will consider a rectangular pulse profile. 

Let T be the time delay between an off-energy particle and the synchronous particle 
which arrives at the center of the bucket t — 0. The equations of motion are 



Here 7] is the slip factor, and Tq the revolution period. Since the bucket is centered at 



di 



AE 




'0 



(1) 
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r = 0, there is no change of energy for — T2/2 < r < T2/2. The second order equation 
of motion for r is 



V 



-eV{T) = 



(3^EoTo 

From the first order equations of motion it follows that the Hamiltonian is 



H 



1 1] 2 6 

'2^^ ' "To Jo 



V{s)ds 



-(AE)^ - -f/(r) 



(2) 



(3) 



where U is the potential function. 

At the bucket boundary or separatrix, the bucket height in energy AEhucket is de- 
termined by 



Hhucket 

Hence the bucket height is 



2p^-C/0 io 



T2/2+T1 



32/2 



bucket 



T2/2+T1 



V JT2/2 



eV{s)ds 



V{s)ds 



1/2 



(4) 



(5) 



On any other level curve inside the boundary, the maximum time extent of a particle 
is T2/2 + W, where the parameter W is the barrier penetration depth for that particle. 
Then the maximum energy extent of a particle as a function of W is 

1/2 



AE 



'2(3^Eo 




rT2/2+W 


To 


V 


JT2/2 



eV{s)ds 



(6) 



The particle energy deviation stays constant at this value for — T2/2 < r < T2/2. At 
other times, the energy deviation can be found from 



AE{t) 



{AEf - 



eV{s)ds 



1/2 



1^71^0 



1 1/2 



eV{s)ds 



T2/2+T 



(7) 



From the first of Equations ([T]) it follows that the time period on a phase curve is 



[2^1 + 4 / 
AE Jo 



dr 



AE{t) 



where AE in the integral is determined from the expression in Equation ([7]). 
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We now specialize to tiie case of rectangular pulse profiles, i.e. the voltage is given 



by 



V{t) = -Vo 
= 
= Vo 



-{It2 + T,)<t<~T2 

2 ^ - - 2 ^ 

-Ta < T < -T2 + Ti 
2 ~ 2 



(9) 



The potential which is only determined up to an additive constant can be written as 



U{t) = -Vo[t+-T2]<0 



= 
= Voir 



zTo] > 



- -Ts < r < -T2 
2 - 2 

It2<t<^T2+T^ 



(10) 



Here we have chosen the additive constant so that the potential also vanishes where 
the voltage does. The potential is zero outside this range defined above. 

For this profile, the peak energy deviation on a phase curve with barrier penetration 
depth W, the bucket height and the time period on an orbit are 



AE 



2f^eVoW 
. \v\To 



nl/2 



bucket 



'\v\To 



eVoT, 



1/2 



T, = 2 



To 



^\ + A^AE 
m AE eVo 



(11) 



Note that the period does not go to infinity on the separatrix, as is the case for harmonic 
rf. 

The peak energy offset, at which the period reaches a minimum, is 



A£;| 



minTs 



P^EoeVoT^ 



1/2 



4:Tl 



-AE, 



bucket 



(12) 



2 1^1 To_ 

This leads to the conclusions that (a) if T2 > 4Ti, the period minimum occurs outside 
the bucket, i.e. the synchrotron period just decreases monotonically inside the bucket 
and (b) if T2 < 4Ti, there is a point where the period reaches a minimum and then 
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increases up to the bucket. The region of zero slope in the period is associated with 
local loss of Landau damping. This simple stability criterion requires that the interval 
between pulses be greater than four times the pulse width. 

The synchrotron tune on an orbit and the maximum tune when it occurs inside the 
bucket are respectively 

-1 



To 



+ 4- 



AE 



\v\ AETq eVo 



T2 + 4W 



\ri\ eK) To 



1/2 



(13) 



.32/?2£;oT2. 

Note that this maximum tune is independent of the pulse width T^. 

The peak energy AE on a phase curve is related to the barrier penetration depth 
W by 

A P, fW 

(14) 



AE 

The synchrotron tune on a phase curve is related to the maximum tune as 

ff^ AE 




K 



K 





(15) 



^s,max 1 + 4aC^ ' V ^2 V ^2 ^E^y^^ket 

where we have defined a dimensionless parameter k to label a phase curve. 

On the phase curve where the synchrotron tune is maximum, k — 1/2. At the 
bucket boundary, it has its maximum value Kbucket = \/^T7^- 

Area under phase curve 

The phase curve for a rectangular pulse profile consists of straight lines at ±AE 
between — T2/2 and T2/2 closed by parabolic paths. On these parabolic paths, the 
relation between the time delay coordinate and the energy deviation is 



2P''EoeVo 

while the area under each parabolic path is 



[{AEy-{AEy], 



(16) 



(17) 
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Hence the area under a phase curve is 

= '^^^^ + lw^:^^'>' = + r^i 

while the bucket area is 

4Ti 

Abucket = '2A.EbucketT2[l + oT^t] = 2 

If the maximum energy spread of the bunch AEbunch is known, e.g. from a longitudinal 
Schottky measurement, then the bunch area can be found from Equation (fT8|) . Cur- 
rently the energy spread and the bunch area of the Recycler beam are measured using 
a 1.7GHz Schottky detector. For a beam with energy spread less than ±23 MeV, the 
Schottky frequency spread is less than ±40 kHz, Hence in practice, the area under the 
Schottky spectrum between -40 MHz and 40 MHz is calculated, and the energy spread 
corresponding to 95% of this area is found. This value is used to calculate the bunch 
area which contains 95% of the particles. 

Bunch length 

The bunch barrier penetration depth Wh can also be determined from the maximum 
energy spread in the bunch using Equation dHl), i.e. Wb = {AEbunch/ AEbucketYTi. The 
bunch length is u = T2/2 + Wb- These equations express the fact that larger initial 
energy spread in the bunch lead to deeper barrier penetration. 



4Ti 



(19) 



III. STATIONARY DISTRIBUTION IN A BARRIER BUCKET 

The time evolution of the phase space density p follows from 

If the density is not explicitly time dependent and it is a function of the Hamiltonian 
H so that the Poisson bracket {p, H} vanishes, then the density is stationary. 

Due to the symmetry of the barrier, the phase curves are also symmetrical, each 
phase curve extending from —W — O.5T2 to ly ± O.5T2 on the r axis (this defines W) 
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and from —AE to on the energy axis. The peak energy on a phase curve can 
be found from 

(A£)^ = ^ V(.)d. = (ir, + W) - £/(1t,)1 (21) 

1^1 -To J0.5T2 IVlTo 2 2 

At the bunch boundary on the r axis t — Tb and AE — 0, and the Hamiltonian takes 

on the value 

H = Hb^ -^U{n) ^ -^Ub (22) 
-to -to 

written in terms of the potential function 1/^ si t — t^. 



A. Binomial distribution 

Consider a general binomial distribution for the density 

p{AE,r)^c,[Hb-H]P (23) 

where p is a real number and Ci is the normalization constant. 
The density can be expanded to the form 

p{AE,r) ^ c,[-^n^^[U{r)-Ub]-AEr 

^ <^i[^mMrr - ^E'r (24) 

where we have used the fact that at any point on the bunch boundary, the energy 
deviation can be found from 

AE,{rr^^-^[U{r)-Ub] (25) 

The line density is obtained by projecting the phase space density on to the r axis, 

A(r) = / p{AE, r)d{AE) = 2c,[-^Y / [AE,{rf - AE^f (26) 

where in the second equality we have used the fact that the integrand is even. Using 
the integration result 



[AE^ - AE^]Pd{AE) = \ ,iA AEb{T)]^P+^ 



I ' ' '"'^ ' r(p + 3/2) 

9 



we obtain for the line density 

A(r) = Xo[U{t) - U,f+'/', Ao = N,[ T {U{t) - U,y+'/']-' (27) 

where we have absorbed the constant Ci and the other parameters into a new constant 
Ao. 

The distortion of the potential due to external inductive impedances and space 
charge can be included with the concept of an effective impedance. The effective 
impedance at frequency cu is defined as 



t[uJoL - — — 

n Wl 



i{ujQL - —5-^) = icooLeff (28) 



where n = uj/uq is the harmonic of the revolution frequency, L is the inductive 
impedance of the machine, Zq = 377fl is the vacuum impedance. The geometric 
space charge factor go for a round beam in a round beam pipe is given by 

go = ^ + 2ln{b/a) (29) 

where b is the radius of the beam-pipe and a is the beam radius. 
The induced voltage due to this impedance is 

Prom the bunch current /(r) = eA(T), it follows that for a binomial distribution, 

f = eAo(p + ^)[C/(r)-C/,r^y(r) (31) 

Introducing the notation 

u{-n, n) = r [u{t) - u,Y+-^dT (32) 

J -Th 



the induced voltage can be written as 



V^na = - '^f +^^ m[^][^(r) - U,r'^V{r) (33) 
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FIG. 1: Rectangular barrier bucket voltage profile and the line density (in blue) in a stationary 
elliptic distribution. Ti is the width of each pulse, T2 is the duration between the pulses, Tb 
is half the bunch length and Wb is the barrier penetration depth. 

B. Elliptic distribution 

The elliptic distribution is obtained by choosing p = 1/2 which yields 



A(r) = Ao[f/(r) - Ub] 



(34) 



As we will see, this distribution has some special properties, first pointed out by Hof- 
mann and Pedersen 7|]. 

The line density for an elliptic distribution with boundaries at {—Tb, Tb) is 

Nb 



A(r) = Ao[f/(r) - f/(r,)] 



-[UiT)-UiTb)] 



U{-Tb, Tb) 

Hence the stationary elliptic distribution in this rectangular barrier bucket is 



A(r) 



[r'b - (1^2)^] 



+ Tb 



Nb 



[r - n 



Tb<T < T2 

-T2 < r < -T2 
2 ~ 2 



-T2<T< Tb 



(35) 



(36) 



This has a trapezoidal shape as shown in Figure [TJ 
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For an elliptic distribution the induced voltage is 

V^n, = f^\m[^]V{r) (37) 

In this case it has the same form as the rf voltage so the total voltage is related to the 
focusing rf voltage by a change of scale. The total voltage is 

V,{r) = V{t) + VUr) = |l ^^^Im[^] W(r) (38) 

By definition of the voltage waveform, a synchronous particle sees no focusing volt- 
age. Hence the ratio of the total focusing voltage to the rf focusing voltage is 

kt = — = l -Im^ (39) 

V LJoU[-Tb,Tb) n 

The maximum synchrotron tune in the presence of space charge and external 
impedances is 

max _ r l"^! ^0il/2 _ . /F max / Ar.\ 

where ly^^^ is the maximum tune in the absence of space charge and external 
impedances. Similarly the area under a phase curve with turning points {ti,T2) on 
the T axis is given by 

A(ri,r2) = v^A(ri,r2) (41) 



The area of the bucket is reduced by the factor ^/^^"cfeet gjygj^ Equation fl5^ with 
Tfo = O.5T2 + Ti. The maximum theoretical intensity that can be stored in the bucket 
corresponds to the case where the space charge and external impedances distort the 
potential sufficiently to reduce the bucket area to zero. In practice, the limiting bunch 
intensity will be lower than this value, for example when the reduced bucket area can 
just accomodate the longitudinal emittance of bunches from the Booster. 



1. Energy distribution 

The energy distribution can be found by projecting the phase space density on to 
the energy axis. Let fi{AE) represent the energy distribution, then 

/x(AE) = J p{AE,r)dT = c^[^^f' £[\AE',{t)-AEY^' (42) 
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Note that /i is even in the energy. The energy offset on the bunch boundary is defined 
in Equation (!25l) . Over the three regions we have 



bucket 



'^■'^ bucket 



Ti 

n - 



^^bucket 



r-n 

Ti 



-n<r< -^T2 

- -Ts < r < -T2 
2 - 2 

-T2 < r < n 
2 - 



(43) 



Consequently 



/i(AE) = ci[-^ 



eK)ii/2 



'-^^bucket 



'-^^bucket 



+ 



T"6 



^-^^ |l/2j 
T2/2 ^-^bucket 



^^bucket Ti 



1/2 



3 



AE 



AEi 



bucket 



^^bucket 



3/2 



(44) 



This has its maximum at AE = and vanishes at 

1 



AE„ 



±AE, 



bucket 



b + 



2a{2 - 3ab) 



[3a%^ - 1 + ^/{3a%^ - If + Aa^h^{2 - 3ah)] 



1/2 
(45) 

where a = (4/3)Ti/T2, b = Wb/Ti. For the Recycler parameters shown in Table [U 
this yields a vanishing of the energy distribution at ±7.87 MeV well below the bucket 
height AEhucket = 17.38 MeV. 

A sketch of the density distribution is shown in Figure [21 the parameters are those 
of the Recycler. The plotted density is normalized to its maximum value at AE = 0. 
As is the case for the line density, there is a well defined value beyond which the 
distribution fi{E) vanishes. 



C. Exponential distribution 



The energy distribution for the elliptic phase space distribution does not match 
the approximately Gaussian energy distribution that is observed in the longitudinal 
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-2 2 
AE [MeV] 



FIG. 2: The energy distribution as a function of the energy spread, for the Recycler param- 
eters, assuming an ehiptic phase space distribution. 

Schottky spectrum of Recycler bunches. We therefore consider another distribution 
which naturally leads to a Gaussian energy distribution. Consider the density to be an 
exponential function of the Hamiltonian. 



p{H) = poexp[H/Ho] = poexp{ — [ 



1 r -V . eU{r) 



]} 



Ho'2l3^Eo' ' To 
The line density obtained by projecting the density onto the time axis is 

ef/(r) 



(46) 



A(r) = y p{AE,T)dAE = Aoexp[ 



HqTq 



(47) 



where Aq is a normalization constant and Hq is a scale constant to be defined later. 
From the expression in Eq (fTOj) for the potential function in a rectangular barrier 
bucket, it follows that the line density for this distribution is 



A(r) = Aoexp[^(r + lr, + rO 

= Aoexp[^r;-7;^TiJ 



Ao exp 



HqTq 



- (Ti + ^Ta) < r < -^Ts 

- -T2 < r < -T2 
2 ~ 2 



,r--T2-Ti)l -T2 < r < -T2 -f Ti 
ifoTo 2 2 2 



(4^ 



Figure [3] shows a sketch of the distribution. 



14 




FIG. 3: Sketch of the hne density for an exponential distribution. Parameters are arbitrary. 

Note that the density does not exactly vanish at r^, instead A(— r^,) = A(r{,) = 
Xoexp[-{eVo[Ti - Wb])/{HoTo)] where as before Wb = n - \T2. If eV^lTi - 
Wb\/ {HqTq) ^ 1, then the density is small and can be effectively taken to vanish 
at r = Tft. 

Since the density does not exactly vanish at ir^,, it would be more accurate to 
extend the region of non-zero density to the entire extent of the barrier bucket. In that 



case the normalization condition is / 



A(r)(ir = A^^^ which yields 



exp[-^^^J) + r2 



(49) 



eVo ' HoTo 

If the line density is known in the flat region, e.g. at the center r = 0, then from 
the expression for the line density in Eq fHS]) above, it follows that the normalization 
constant is known, i.e. Aq = A(r = 0) and the normalization condition can be used 
to find the scale constant Hq. It is simpler here to determine this constant from the 
energy distribution which in this case is a Gaussian 

(AE)2 



/i(AE) = / dTp{AE,T) = /ioexp[- 



2a| 



(50) 



where the rms energy deviation a"^ = P'^EqHq/i]. If the rms energy deviation is mea- 
sured from the Schottky spectrum, then Hq is known. 
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The line density written above is not self-consistent as it does not take into account 
the potential well distortion due to space charge. Including this intensity dependent 
effect, the line density is 

A(r) = Aoexph^^] (51) 

Ut{T) = Urf{T)+ r V,^a{s)ds = Urfir) + —lm[^][X{0) - X{t)] 
Jo ^0 n 

It follows that the density at time r = is A(0) = Aoexp[— ef/r/(0)/(i^o^o)]- Since we 
have defined the rf potential such that Urf{0) = 0, it follows that A(0) = Aq. Hence 
the self-consistent density is the solution of the following equation 

A(r) exp[--^Im[^]A(r)] = Ao exp[-^Im[^]Ao] exp[-^^] (52) 

This is the equivalent of the Hassinski equation which defines a self-consistent station- 
ary solution for a sinusoidal rf and an arbitrary wake field. This equation can be solved 
numerically to find the self-consistent density at a given intensity. At low intensities 
this solution will reduce to the form in Equation fj48|) . 



IV. FREQUENCY OF COHERENT MOTION 

Imagine that the bunch is displaced from its center by Ar as a result of which it 
starts to perform small amplitude oscillations. If we let A(r) be the line density of the 
displaced bunch, then it follows that A(r) = A(r — Ar). The infinitesimal force on a 
slice of thickness 5t is 

dF = -^\{t - Ar)dT^ (53) 
pc or 



The total force on the bunch is 

e 



F = -j- j X{r-AT)—dT 



= ^Ar j X'{T)—dr + 0[{Arr] (54) 

where we have expanded the displaced density to first order in Ar and used the fact 
that the average force on the undisplaced bunch vanishes. From the equation of motion 
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and the force law F = — e\^(r)/(/5c) = mefff, it follows that the effective mass of a 
particle is w^e// = PEoTo/^cr]) and hence the effective mass of the bunch is m^^yj'^^ = 
Nh(3 EqTq / {ct]) . If the centroid is oscillating at a frequency Uc, then from the small 
amplitude equation of motion F/mlY}'^^ = = it follows that the coherent 

frequency is 

This coherent frequency is independent of intensity. It turns out that for parabolic 
bunches in a single harmonic rf potential, the coherent frequency and the small ampli- 
tude incoherent frequency coincide. In general these frequencies are not the same for 
other distributions and other rf systems. 



A. Landau damping for an elliptic distribution 

We will consider the rigid dipole mode to be Landau damped if the frequency of 
the rigid dipole mode is within the band of incoherent synchrotron frequencies and 
undamped if it lies outside this band. First we evaluate the coherent dipole frequency 
for an elliptic distribution. Using the line density of a matched binomial distribution, 
the coherent frequency using Eq. is 



This frequency, in the absence of external forces, is independent of the bunch intensity. 
For the elliptic distribution this reduces to 

2 e\ri\ n ^ , , 2eVo\v\ 



'^^ - p^E,nu{-n, n) ^ ^^^^^ - p^E,Un + \T,) ^^^^ 

As remarked above, the total focusing voltage is changed by space charge and ex- 
ternal impedances. For an elliptic distribution, the total voltage is related to the rf 
voltage by the factor kt calculated above and the synchrotron frequency is related to 
the bare synchrotron frequency by yfkf In a barrier bucket, the incoherent frequency 
as a function of amplitude rises from zero at the origin to a maximum value which may 
occur at an amplitude within the bucket. The maximum bare incoherent synchrotron 
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angular frequency cuTo*^ "^^^ ^le found from Equation f[T5]) . In the presence of space 
charge and external impedances it is modified to 

, ,max _ rr, max _ Ti , eA^bIm(.Z'e///n) y^j^ ^^.x ( r:o\ 



In terms of the bare maximum frequency w^q^ , the coherent frequency can be written 

To 



as 

4 T 

r ^ 1I/2 max /cnN 

= -[ , 1^ ] ' ^s,o (59) 

Tf, + 2^2 

Writing rf, = T2/2 + Wb, where Wb can be considered as the barrier penetration depth 
of the bunch, the requirement that the coherent frequency be less than the maximum 
bare incoherent tune is Wb/T2 > — 1. 

Below transition with space charge as the dominant impedance, the net focusing 
voltage is reduced and the incoherent frequency decreases with intensity. At the thresh- 
old for the loss of Landau damping, the maximum incoherent frequency equals the co- 
herent dipole frequency and falls below it at higher intensities. The threshold intensity 
for the loss of Landau damping found by equating cj^"^ and ujc is 

^thresh _ ^0'^0(^b ~ (1^2)^^) . A .2 ^2 . , . 

~ eIm(Z.,,/n) V ^^Tl^^ 

Putting in numbers for the Recycler, we find that 

"''^ L25 (61) 



, ,max 
s 



This says that even at zero intensity, the coherent frequency would be outside the 
incoherent band. As the intensity increases, the incoherent frequencies will decrease 
and the ratio above will increase. Hence there would be no Landau damping at any 
intensity. We conclude that this distribution cannot be suitable to describe the Recycler 
longitudinal density. 



18 



B. Landau damping for an exponential distribution 



The coherent frequency can be found using Eq flSSl) except that we now extend the 
region of integration so 

2 e\r]\ Z-^^^+Ti ^^^_^dU 



We find 

'"^ /?^EoTo2i^f (l-exp[-f|])+T2 ^''^ 
Expressed in terms of the maximum of the bare incoherent frequency uj'^q^, 

, x = P'EoeVo/i\r]\alTo) (64) 



_ 4 1 - exp(-xWfe) ,1/2 



O0so,max 71^1 + (2/(xT2))[l - exp(-xl^;,)] 
Putting in the Recycler numbers we find that in this case 

UJr. 



,max 

This distribution is also unsuitable. 



1.27 (65) 



V. MEASURED DISTRIBUTIONS 



During a study in April 2009, protons were injected into the Recycler with a in- 
creasing number of Booster batches. The bunch profiles were measured with a wall 
current monitor at three different intensities. Longitudinal Schottky spectra were also 
recorded. The measured bunch length profiles are well described by the function 

A^,i(r) = A'ba {[1 + tanh(6r + c)]9{-r) + [1 - tanh(6r - c)]9{t)} (66) 
\fa{0) = A^ba(l + tanhc) 

Here 6 is the step function, a, b, c are the three parameters of the fit. a sets the overall 
normalization for the density, 6 is a measure of the slope in the rising and falling 
regions and the width of each half section is directly proportional to c. Note that the 
same parameters describe both the head and the tail of the profile, implying perfect 
symmetry about the center. 
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FIG. 4: Comparison of the measured line density and fit to the data with a tanh profile at 
three different intensities . 
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FIG. 5: Measured Schottky power of a proton beam in the Recycler compared with a Gaussian 
fit. 
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The measured profiles and the fit function are shown for all three sets in Figure ID 
This fitted function appears to describe the measured profiles with very good accuracy. 
Note that the same fit works equally well at all the intensities suggesting that distortion 
due to intensity effects are not yet important at these levels. This also suggests that 
the tanh fit is probably not the result of a self-consistent solution to the "Hassinski 
equation" flS21) . the assumption being that the fit is also good for lower intensities 
than those measured here. There is some bunch lengthening as the intensity increases 
showing the impact of the capacitive space charge impedance. 

The line density can also be calculated as a function of the incoherent synchrotron 
tune. Figure [6] shows the density relative to the density at the center as a function of 
the bare incoherent tune. The Main Injector ramping frequency is also shown in this 
figure; beam particles resonant with this frequency can be driven to larger amplitudes 
and lost. However the relative density at this frequency is about 2%, the net effect due 
to this resonance should be small at present parameters. 

We can now use this form of the line density to calculate the coherent frequency. 
From the expression for the coherent frequency in Equation fl55|) 



= t\ {-\K-\t2) - K-n)] + [Kn) - \{\t,)]] 

P^'EoToNb 2 2 

2 \r]\ (^Voa^^^^^^^^ _ ^) tanh(-6T2 - c)} (67) 



Hence the ratio of the coherent frequency to the maximum of the bare incoherent 
frequency is 



— a/oT^ -i /tanh [fer^ — c] — tanh[-6T2 — c] (6^ 



V ^ I / I" ■ u "J - ^ L r\ 

Here we use the value of the parameter a determined by the normalization condition, 
i.e. 

1 r 1 ^ 1 , / cosh c , , 1 , , 

'■ = 5l5^- + ^'n'°g< cosh|Kr./2 + ro-c] " («^' 

Using the values of the fit parameters b, c for the data measured in the April 2009 study 
we find that 

= 0.93 (70) 
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Density as a function of tune 
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FIG. 6: The relative density as a function of the bare incoherent tune. Most particles have 
zero tune. For comparison, the tune corresponding to the frequency (0.45 Hz) of the Main 
Injector ramp is also shown. 

i.e. at zero intensity the coherent frequency is within the incoherent spread. Figure 
[7] shows the coherent tune for the three distributions considered here and the bare 
incoherent tune as a function of the parameter W. 

A. Threshold for loss of Landau damping 

We can now take into account the potential well distortion due to space charge and 
find the incoherent tune in the presence of space charge. The total voltage for an 
arbitrary distribution is 



Vt = Vo- — 




(71) 
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Bare incoherent tune and coherent tunes 
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FIG. 7: Bare synchrotron tune vs W, the barrier penetration parameter for each orbit. Also 
shown are the coherent tunes of the dipole mode for three distributions. The coherent mode 
would be damped only for the Tanh distribution. 

The energy deviation on a curve with barrier penetration W at time r is 



{AEf{r;W)- 



JT2/2+T 



eVt{s)ds 



I /TT. X e ^ ,Zeff... ,1„ 1, 



I ^ , VoiW - t) Im(^)[A(-T2 + W) -X{-T^ + r)] (72) 



Introducing the function f{W) and the parameters C, d as 

e ^ ,Zf,ff,_ , 1 



f{W) = W + C^a.n\i{hW + d), ( 



UJqVo 



lm{^)Nba, d = -6T2 - c 



n 



we can write 



2f3^eVoEo 
1^71^0 



IfiW) - fir)] 



(73) 



(74) 



while the peak energy is found from = AE{0, W). The intensity dependence is 
contained in the parameter (. 

As the intensity increases, the bucket height and bucket area decrease. A threshold 
for the maximum intensity that can be stored in the bucket is set when the bucket area 



23 



Bucket height and area vs intensity 
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FIG. 8: The bucket area and bucket height as a function of the bunch intensity. Over the 
range of intensities plotted, the bucket values of the area and height stay well above the initial 
bunch area and bunch energy spread. 

falls to the initial bunch area. Figure [8] shows that at intensities even up to 200 times 
larger than present, the bucket area and height are larger than the initial bunch area 
and energy spread respectively. The threshold is at an intensity of 2.3x10^^, where the 
bucket height equals the initial bunch energy spread of 12 MeV. 

The incoherent synchrotron tune as a function of W can then be found using Eq 



(ED as 



\r]\ 



To 



[f{W) - Ctanhc/]i/2 



w 



dr 



-, -1 



[/(VT) -r-Ctanh(6r + rf)]i/2 

(75) 

The integral in the above expression can be performed numerically. For comparison, 
the bare synchrotron tune without space charge is given in Equation f|T3l) . 

The top plot in Figure [9] shows the relative difference of the incoherent tunes with 
and without space charge at intensities varying over two orders of magnitude . The 
relative tune depression due to space charge is the largest at the inner edge of the voltage 
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Relative incoherent tune shift at various intensities 





n 
u 




-U.i 




-U.Z 




-U.i 




-0.4 




-0.5 




-0.6 




-0.7 




-0.8 




-0.9 




0.4 




3 




U.Z 


b 


A 1 

0.1 







> 




1 

X 


-0.1 


6 


-0.2 


> 




-0.3 




-0.4 




-0.5 



Nb 












Tl = 48, T2= 


324 


/ Nb 


= 3.2E14 




1 




bunch area = 


160.5 eV-sec - 
1 1 



5 10 15 20 25 30 35 40 45 50 

W 

Incoherent tune with space charge - coherent tune 



n 1 1 1 1 1 1 1 r 



Tj = 48, T2 = 324 

Bunch area = 160.5 eV-sec 



_l I I I I I I L. 



123456789 10 
Bunch intensity [10^"^] 

FIG. 9: Top: Relative difference of the incoherent tune with space charge and without 
space charge as a function of W, the barrier penetration at bunch intensities varying from 
A^b = 4 X 10^2 ^ 3 2 X 10^^. Bottom: Difference of the maximum incoherent tune 

with space charge and the coherent frequency (multiplied by 10^) as a function of the bunch 
intensity. The zero crossing defines the threshold for the loss of Landau damping, here the 
threshold is iV;, = 4.7 x 10^^. 
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pulse since the tune is also smallest there. The incoherent frequency decreases as the 
intensity increases and at a threshold frequency, the maximum of the incoherent tune 
will will equal the coherent frequency. At higher intensities, the coherent frequency will 
be outside the spread and Landau damping is lost. For the nominal value of T2 = 324 
'bkts', this threshold is reached when the intensity exceeds 4.7x10^^, about 100 times 
higher than present intensities. 

VI. INSTABILITY THRESHOLDS AT DIFFERENT BUNCH LENGTHS 

It has been observed in the SPS that a region where the incoherent frequency as a 
function of amplitude has an extremum is associated with the appearance of a local 
instability jsl. In a barrier bucket, the maximum of the incoherent tune occurs when 
the ratio of the peak energy on the orbit to the bucket height is A/T'2/ (4Ti), see Equa- 
tion f|T2|) . If this ratio is greater than unity, as is the case for the nominal Recycler 
parameters T2/ (4Ti) = 324/ (4 x 48), the maximum lies outside the bucket. The effects 
of this local instability are not present under normal operation but can be studied by 
changing the separation T2. 

The energy acceptance or bucket height depends on the area under the voltage pulse 
VqTi but not on T2. If we change the bunch length by adiabatically changing T2, the 
bunch area will be preserved and the energy acceptance will also be constant. The 
bucket area however is determined by the bucket height and the value of T2. The value 
of T2 when the bucket area just encloses the bunch area is 



Lower bunch area leads to a lower value for T™*" and hence allows a wider range of 
variation with T2 < 4Ti. 

These studies were carried out in the Recycler using protons. The beam cooling 
systems as well as the transverse dampers were turned off during the experiment. For 
these measurements, a bunch with much lower intensity than normal was injected - 
the initial bunch area was about 27 eV-sec compared to a more typical value of 100 




1 



(76) 
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eV-sec. The bunch was injected into a gap which was 84 Booster buckets long, i.e. 
Ta = 84 'bkts' (a 'bkt'= 0.01893/isec). We chose Tl = 48 bkts and pulse height 
of 1.8kV, these being the usual values during regular operation. Once equilibrium 
was reached after injection, beam intensity along with data from the wall current and 
Schottky monitors were recorded to establish the initial parameters. Subsequently the 
beam was expanded or compressed adiabatically by changing T2 to different values of 
interest without changing Ti. 

Figure [TO] shows the measured longitudinal profiles along with the fits by the tank 
function. In all cases, the fit describes the measured profile quite well. One slight 
exception is for the shortest bunch with T2 = 14 bkts where the fit at the center 
does not quite match the flatness of the measured profile. Furthermore, it cannot 
be explained by the wall current monitor saturation which occurs at a minimum of 5 
times the bunch intensity used in this experiment. However the rising and falling edges 
match very well. These edges in the barrier pulses determine the coherent frequency, 
not the bunch profile in the center. 

Assuming that the bunch area is preserved, compressing the bunch length increases 
the energy spread and conversely. Figure [TT] shows the maximum energy spread deduced 
from the Schottky spectrum and the theoretical values (shown by the curve) assuming 
that the initial bunch area was conserved. There are discrepancies between the two, 
especially at the shortest T2 = 14 bkts. Some of this can be attributed to measurement 
errors of the energy spread from the Schottky spectrum. It is also likely that several 
beam dynamics processes lead to an increase in bunch area and loss of particles as the 
gap T2 was decreased. This is confirmed by the measured Schottky spectra for the 
different profiles shown in Figure [121 At the smallest values of T2 , the energy spread 
exceeds the energy acceptance of 17.4 MeV. However the spectral amplitudes, shown 
in db, also show that the losses were of the order of 1%. 

We now consider the coherent frequency and the ratio of this frequency to the 
maximum of the bare incoherent frequency as a function of the bunch length, using the 
energy spread and the fitted values of the parameters b, c for each of the bunch profiles 
shown in Figure fTOl This ratio is shown in Figure fT3l First, the ratio increases with 
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FIG. 10: Measured bunch profiles in March 2009. The value of T2 was changed while Ti was 
kept constant. The sequence of T2 values is the measurement order. The measured data and 
the fits with the tanh function are shown for the eight profiles. 
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Maximum energy spread of bunch 




T2 [bkts] 

FIG. 11: The maximum energy spread as a function of T2. The data points are the values 
obtained from the Schottky spectrum while the curve shows the expected value assuming 
that the initial bunch area is preserved and no beam is lost. Over this range of T2 from 
14 bkts to 324 bkts, the maximum spread stays below the energy acceptance of 17.3 MeV. 
The deviations between the data and the curve show that the bunch area was not preserved, 
especially at the smallest value of T2 = 14bkts. 

the bunch length or T2 before falling at the largest value of T2. Note that both the 
coherent frequency and the maximum of the bare incoherent frequency decrease with 
T2 ■ oJsQ^max decreases as 1/ a/T^ while the coherent frequency decreases at a slower 
rate, hence the ratio increases with T2. One would conclude from this observation that 
for bunches with the same (or nearly the same here) longitudinal emittance, shorter 
bunches have a higher threshold for the loss of Landau damping. Secondly, we observe 
that except for T2 = 14, this ratio is greater than one implying that Landau damping 
was lost for most of the bunches in this study. This may explain some of the observed 
increase in bunch area. 

Fig. [14] shows the surviving beam as a function of time. Various T2 to Ti ratios 
along with the beam lifetime are also indicated. There are three distinct regions in 
this plot. The first section had T2 < 4Ti with T2 assuming values of 14, 44, 84, 130, 180 
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FIG. 12: Schottky spectra of bunches with different values of T2 after background subtraction. 
At the smahest values of the spectra extends beyond the energy acceptance suggesting 
that there was some beam loss but this amount was of the order of 1%. 
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FIG. 13: Ratio of coherent to maximum of the bare incoherent frequency for the different 
bunches with line profiles shown in Figure [TOl 
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bkts. The average lifetime of the beam at 8 hours was significantly lower than in the 
second section where it was about 50 hours with T2 = 228, 276 and > 4Ti in both cases. 
In the third region the beam was compressed to its original value of T2 = 84 bkts. The 
lifetime again decreased by about a factor of two. There are many possible sources for 
these changes in lifetime as T2 was changed: different energy spread, loss of Landau 
damping, local instability when T2 < 4Ti, and possible resonance with the frequency of 
the Main Injector ramping. The density of particles at this resonant frequency changes 
with the bunch length and this resonance has been known to cause beam loss in the 
past. While we do not at present know the relative importance of these effects, the 
data is not inconsistent with the hypothesis that a local instability at the extremum 
of the incoherent synchrotron tune can contribute to beam loss. The effects of this 
instability are easily avoided in the barrier bucket by choosing appropriate values of 
Ti,T2 




45 

Time (min) 



FIG. 14: Lifetime at various times with different T2 to Ti ratios. 



31 



VII. MICROWAVE INSTABILITY 



So far we've only discussed the effect of the capacitive space charge impedance. 
However the Recycler also has a broadband resistive impedance with a cutoff at /c ~ 
c/h 3.95 GHz. At sufficiently high intensities, this resistive impedance can lead to 
asymmetric bunch shape distortion and eventually to turbulent bunch lengthening due 
to the onset of the microwave instability. While there is as yet no satisfactory theory 
of this instability for bunched beams, estimates for the threshold of the instability are 
made by using the coasting beam theory and replacing the beam current by the peak 
bunch current. For Gaussian bunches, the intensity threshold is given by 
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We observe that for constant bunch area A ~ (JzCTs, the threshold A^^ oc A^/az- Hence 
for bunches of the same area, shorter bunches have a higher threshold. 

As a preliminary estimate for the threshold of this instability, we will use the rms 
values of the energy spread and the bunch length for the measured distribution during 
normal operation. The rms energy spread is found from the Schottky spectrum while 
the rms bunch length can be found from 

Trms = Vo-l (1 + tanh(6r + c) )(ir + / t^(1 — tanh(6r — c))(i 

i-(T2/2+Ti) Jo 

For the fitted values of the parameters, we find r^ms = 3.095;U seconds while the 
rms energy spread is a^; = 4MeV. The effective wall impedance is estimated to be 
Z\\/n = 1.6 Ohms. Substituting the values we find that the bunch intensity threshold 
for the microwave instability is 

Nf, = 2.95 X 10^^ (78) 
This is about three orders of magnitude above present intensities. Space charge below 



transition does not induce the microwave instability but some results, e.g. {g] suggest 
that it may moderate the growth of the instability. Simulations may be necessary to 
understand this process. 
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CONCLUSIONS 



ur primary findings are as follows: 

We considered two typical stationary phase space distributions as possible candi- 
dates for the Recycler: an elliptic distribution and an exponential (in the Hamil- 
tonian) distribution. The first has an energy distribution which does not match 
observations while the second has a longitudinal profile distribution which does 
not match observations. In addition, both of these distributions are above the 
threshold for loss of Landau damping at present parameters. 

The longitudinal distribution that describes the Recycler bunches is a tanh dis- 
tribution. This is found to be in very good agreement with observations made on 
separate days and with different beam parameters. This distribution was used 
to find the intensity at which the coherent dipole frequency is at the edge of the 
incoherent frequency distribution in the presence of space charge. This threshold 
for the loss of Landau damping is about two orders of magnitude above present 
intensities. 

An experiment to determine the impact of a local instability possibly induced at 
shorter bunch lengths was discussed. There are at least two competing mech- 
anisms operating in this regime. On the one hand, at these bunch lengths the 
extremum of the incoherent frequency lies within the bunch frequency distribu- 
tion and a local instability may develop. On the other hand we found that the 
threshold for the loss of Landau damping increases as the bunch is shortened 
with the area preserved. Nevertheless, beam loss was observed at shorter bunch 
lengths. Possible causes include higher energy spread, resonance with the Main 
Injector ramp frequency and perhaps the local instability. 

A rough estimate of the microwave instability threshold shows that in the absence 
of space charge, the threshold is about three orders of magnitude higher than 
present intensities. Table |TT] shows the intensity thresholds from different effects. 
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Effect 


Intensity threshold 


Loss of Landau damping 
Bucket height shrinks to 12 MeV 
Microwave instabihty 


4.7x10^^ 
2.3x10^5 
2.9x10^5 



TABLE IL Intensity thresholds for different effects, assuming Ti = 48 bkts, T2 = 324 bkts. 

In a subsequent paper we will examine the stability diagrams in the presence of 
space charge and a resistive impedance. Wc will also report on results of simulations 
with a macro-particle tracking code (ESME) and a Vlasov code. The Vlasov code will 
be employed to compare with the intensity thresholds found in this paper as well as to 
understand the bunch behaviour following the onset of instability. 
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